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[l]  On  the  night  of  December  3rd,  2004  (UT  day  338),  we 
observed  a  significant  acceleration  of  horizontal  wind  near 
100  km  between  0900  and  0915  UT  accompanied  by  a 
temperature  cooling  at  the  same  altitude  and  warming  below 
it.  The  Lomb  spectrum  analysis  of  the  raw  dataset  revealed 
that  a  gravity  wave  with  1 . 5  hr  period  was  significant  between 
0500  and  0900  UT,  but  blurred  after  0900  UT,  suggesting  the 
transfer  of  wave  energy  and  momentum  from  wave  field  to 
mean  flow.  Most  likely,  this  observed  phenomenon  is  due  to 
the  breaking  of  an  upward  propagating  gravity  wave  with  an 
apparent  period  of  ^  1 .5  hr.  Using  linear  saturation  theory  and 
assuming  a  monochromatic  wave  packet,  we  estimated  the 
characteristics  of  breaking  gravity  wave,  eddy  diffusion 
coefficient,  and  a  simple  relation  between  Prandtl  number 
and  turbulence  localization  measure  when  the  wave  is 
breaking,  from  the  experimentally  determined  heating  rate, 
horizontal  wind  acceleration,  and  background  temperature 
and  winds.  Citation:  Li.  T.,  C.-Y.  She,  H.-L.  Liu,  and  M.  T. 
Montgomery  (2007),  Evidence  of  a  gravity  wave  breaking  event 
and  the  estimation  of  the  wave  characteristics  from  sodium  lidar 
observation  over  Fort  Collins,  CO  (41°N,  105°W),  Geophvs.  Res. 
Lett.,  34,  L05815,  doi:10.1029/2006GL028988. 

1.  Introduction 

[2]  Atmospheric  gravity  waves,  since  first  observed  and 
studied  in  the  ionosphere  by  Hines  [1960],  have  been  a 
major  research  topic  for  many  decades  because  of  their 
dramatic  impacts  on  the  energy  and  momentum  budgets  in 
the  middle  and  upper  atmosphere  [Fritts  and  Alexander, 
2003].  Because  the  wave  energy  density  is  constant  for  a 
gravity  wave  propagating  energy  upward,  the  amplitude  of 
the  wave  grows  with  altitude  to  compensate  the  exponential 
decease  of  air  density.  In  the  mesopause  region  (80-105  km), 
the  gravity  waves  could  become  unstable  and  tend  to  break 
via  either  convective  or  dynamic  instability,  and  thus  turbu¬ 
lence  is  generated  within  the  wave  field  [Fritts  and  Rastogi, 
1985],  The  breaking  gravity  wave  will  deposit  momentum 
into  the  mean  flow  and  induce  vertical  heat  flux,  leading  to 
wind  acceleration/deceleration  and  temperature  changes 
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around  the  breaking  region  [Lindzen  1981;  Liu  and  Hagan, 
1998], 

[3]  In  recent  years,  many  gravity  wave  breaking  cases 
were  studied  with  the  temperature  vertical  profiles  observed 
by  Rayleigh  lidar  [Hauchecorne  et  al.,  1987;  Meriwether  et 
al.,  1998]  and  sodium  lidar  [Williams  et  al.,  2002].  But  there 
are  very  few  publications  on  the  gravity  wave  breaking 
revealed  by  simultaneous  observations  of  both  temperature 
and  horizontal  wind  with  high  vertical  and  temporal  reso¬ 
lutions.  The  Colorado  State  University  (CSU)  sodium  lidar 
system  located  in  Fort  Collins,  CO  (40.6°N,  105°W)  is  one 
of  the  most  robust  ground-based  lidar  systems  and  is 
capable  of  simultaneously  observing  temperature  and  hor¬ 
izontal  wind  in  the  mesopause  region  on  a  24-hour  contin¬ 
uous  basis,  weather  permitting.  Many  valuable  datasets 
have  been  collected  under  this  operational  mode  since 
May  2002  [She  et  ah,  2004],  In  this  paper,  we  present  lidar 
data  from  the  night  of  December  3rd,  2004,  and  the 
observation  of  a  large  sudden  horizontal  wind  acceleration 
near  100  km  accompanied  by  a  temperature  change  impli¬ 
cating  a  sudden  downward  heat  flux  transport.  These 
observed  features  are  highly  suggestive  of  wave  breaking 
in  these  altitudes.  With  the  observed  heating  rate  and  wind 
acceleration,  we  estimated  the  characteristics  of  the  break¬ 
ing  gravity  wave,  using  the  linear  saturation  theory.  The 
estimated  parameters  are  consistent  with  an  observed  wave 
with  a  period  of  ~  1 .5  hr,  whose  spectrum  power  was  greatly 
reduced  after  the  sudden  horizontal  wind  acceleration  and 
temperature  cooling  near  ^100  km  and  warming  below  it. 
Discussions  will  be  presented  followed  by  a  conclusion. 

2.  Sodium  Lidar  Observations 

[4]  The  CSU  sodium  lidar  system  currently  utilizes  a  three- 
frequency  design  that  enables  simultaneous  observations  of 
temperature,  zonal  wind,  meridional  wind  and  sodium  density 
in  the  mesopause  region  on  a  24-hour  continuous  basis.  The 
system  receiver  employs  two  14-inch  Celestron  telescopes 
pointing  eastward  and  northward,  both  at  30°  from  zenith.  The 
lidar  transmitter  directs  two  laser  beams  into  the  atmosphere, 
each  aligned  parallel  to  one  telescope,  with  north  and  east 
beams  separated  by  ~70  km  in  the  mesopause  region.  In 
comparison  to  other  lidar  systems,  our  sodium  lidar  system  has 
a  modest  power-aperture  product  of  0.05  Wm2  per  beam. 
However,  due  to  the  resonance  enhancement  in  laser  induced 
fluorescence,  the  received  signal  is  sufficient  to  probe  atmo¬ 
spheric  instabilities  as  well  as  various  atmospheric  waves  and 
their  interactions  [She  et  ah,  2004;  Li  et  ah,  2005],  The 
precision  of  the  measurements  for  temperature  and  wind  with 
a  2  km  vertical  resolution  and  15  min  temporal  resolution 
under  the  nighttime  clear  sky  condition  is  estimated  to  be 
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(a)  Vertical  profiles  of  temperature,  UT  day  338,  2004 


(c)  Vertical  profiles  of  meridional  wind,  UT  day  338,  2004 


Figure  1.  Vertical  profiles  of  (a)  temperature,  (b)  zonal 
wind,  and  (c)  meridional  wind  between  0800  and  1300  UT 
on  the  night  of  December  3rd  (UT  day  338),  2004.  These 
profiles  are  spaced  at  15  min  intervals;  whose  horizontal 
scale  also  correspond  to  20  K  in  temperature  or  40  m/s  in 
wind  speed.  The  vertical  dotted  lines  denote  200  K  for 
temperature  and  0  m/s  for  wind. 

1.5  K  and  3  m/s,  respectively,  at  92  km  (sodium  peak),  and 
6  K  and  12  m/s,  respectively,  at  both  81  km  and  105  km 
(sodium  layer  edges). 

[5]  During  our  December  2004  campaign,  we  observed 
significant  horizontal  wind  acceleration  and  temperature 
cooling  near  100  km  and  temperature  warming  near  96  km 
between  0900-0915  UT  on  the  night  of  December  3rd, 
2004.  The  vertical  profiles  of  temperature  (a),  zonal  wind 
(b),  and  meridional  wind  (c),  with  15  min  temporal  resolu¬ 
tion  and  2  km  vertical  resolution  between  0800  and  1300  UT  | 
for  this  night  are  shown  in  Figure  1.  These  15-min-averaged  -a 
profiles  are  spaced  at  1 5  min  interval,  which  also  correspond  ^ 
to  20  K  for  temperature  and  40  m/s  for  wind.  The  vertical 
dotted  lines  represent  200  K  for  temperature  and  0  m/s  for 
wind.  On  the  temperature  plots,  we  could  clearly  see  that  the 
temperature  suddenly  warms  up  near  96  km  and  cools  down 
near  100  km  at  around  0900  and  0915  UT.  After  0915  UT, 
this  thermal  pattern  of  cooling  and  heating  started  to  progress 
downward  with  a  speed  of  ~1  km/hr.  In  the  zonal  and 
meridional  wind  components,  we  saw  horizontal  wind 
acceleration  near  100  km  toward  northeast  direction  at 
0900  UT.  After  0900  UT,  zonal  wind  remains  large  in  the 
eastward  direction  at  1 00  km  for  several  hours  without  clear 


downward  progression.  The  positive  peak  of  meridional 
wind  at  ~97  km  grew  much  larger  after  0900  UT  than  the 
peak  at  the  same  altitude  before  0900  UT  and  progressed 
downward  slowly.  However,  the  negative  peak  below  90  km 
progressed  downward  with  a  faster  speed  before  and  after 
0900  UT.  All  these  indicate  the  complexity  of  wave  dynam¬ 
ics  during  this  night. 

[6]  Shown  in  the  Figure  2  are  the  plot  of  zonal  wind  vs. 
UT  hour  at  100  km  (left  panel)  and  the  contour  plots  of 
Lomb  power  spectra  between  0500  and  0900  UT  (right  top), 
between  0900  and  1300  UT  (right  bottom)  based  on  the 
same  dataset  shown  in  Figure  1.  The  90%  significant  level 
is  also  plotted  in  Figure  2  with  the  white  solid  lines.  As 
clearly  shown  in  the  left  panel,  the  wave  oscillation  with 
~1.5  hr  period  and  ~15  m/s  amplitude  persists  between 
0500  and  0900  UT,  which  is  also  present  in  the  Lomb  power 
spectrum  (right  top  panel).  After  0900  UT,  the  mean  zonal 
wind  increases  dramatically,  and  at  the  same  time  the  1.5  hr 
wave  oscillation  disappears.  The  disappearance  of  this  wave 
is  also  confirmed  in  the  Lomb  power  between  0900  and 
1300  UT  at  right  bottom  panel.  Therefore,  it  is  most  likely 
that  the  temperature  warming  and  mean  flow  acceleration 
between  0900  and  0915  UT  are  due  to  the  breaking  of  an 
upward  propagating  gravity  wave  with  an  apparent  period 
of  ~1.5  hr.  Lomb  power  spectra  of  both  temperature  and 
meridional  wind  were  similarly  examined,  the  presence  and 
absence  of  1.5  hr  wave  before  and  after  0900  UT  is  also 
evident. 

3.  Results  and  Discussion 

[7]  The  Lindzen’s  linear  saturation  theory  has  been  the 
base  of  most  GW  parameterization  schemes  used  in  both 
mechanistic  and  General  Circulation  global  scale  Models 
[e.g.  Holton,  1982;  Alexander  and  Dunkerton,  1999;  Sassi 
et  al,  2002],  and  the  applications  have  been  quite  success¬ 
ful.  It  has  been  compared  with  mesoscale  nonlinear  numer¬ 
ical  results,  and  reasonable  agreements  were  found  [Liu  et 
al.,  1999].  Therefore,  we  are  going  to  use  this  theory  to 
estimate  the  characteristics  of  breaking  gravity  wave. 
According  to  linear  saturation  theory,  the  mean  flow  accel- 


Figure  2.  (left)  The  plot  for  zonal  wind  vs.  UT  hour  at 
1 00  km  and  the  contour  plots  of  Lomb  power  spectra  (top 
right)  between  0500  and  0900  UT  and  (bottom  right)  between 
0900  and  1300  UT.  The  90%  significant  level  is  also  plotted 
with  white  lines  in  the  Lomb  power  contours. 
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eration  in  the  region  where  the  monochromatic  gravity  wave 
breaks  is  given  by  [Holton,  1982] 


du 

Q-j-  '  Z=Z break 


=  -jMPoU'w') 


N2Km 

[u-c) 


(1) 


where  p0  is  the  air  density;  u  is  horizontal  mean  wind 
projected  in  the  wave  propagation  direction;  u'w'  is  the 
vertical  flux  of  horizontal  momentum  with  u'  and  W  being 
perturbations  of  horizontal  wind  and  vertical  wind, 
respectively;  c  is  horizontal  phase  speed  of  the  gravity 
wave  (with  c -u  the  intrinsic  wave  phase  speed);  N  is  Brunt- 
Vaisalla  frequency  of  the  background  flow  (including  tidal 
components);  z  is  altitude;  zbreak  is  the  altitude  where  wave 
breaks;  and  Km  denotes  the  eddy  diffusion  coefficient, 
which  is  expressed  as  [Lindzen,  1981]: 


^  k(u  —  c)4  /  1  3  1  du\ 

m  ~  N*  ~~  2u  —  cdz ) 


(2) 
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Figure  3.  The  heating  rate,  zonal  wind  acceleration,  and 
meridional  wind  acceleration  at  0900  UT  on  the  night  of 
December  3rd  (UT  day  338),  2004.  The  breaking  region  is 
near  100  km. 


where  k  is  the  horizontal  wave  number;  H  is  the  scale 
height;  and  du/dz  is  the  vertical  shear  of  horizontal  mean 
wind  in  the  wave  propagation  direction.  Since  the  gravity 
wave  perturbation  is  relative  to  the  background  atmosphere 
at  break  region,  there  is  no  explicit  wave  amplitude 
information  in  the  equation  (1).  The  maximum  heating  rate 
below  the  breaking  region  is  given  by  [Liu,  2000] 


max  |z  <  Zbreak 


To|z 


gl-M  1 


W 

4H 


with  £ 


(3  -  4r;) 
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where  T0  is  mean  temperature  over  the  breaking  region;  r/  is 
a  measure  of  the  localization  of  turbulence;  g  is  the 
gravitational  acceleration;  Az  is  the  vertical  wavelength  of 
the  gravity  wave;  Pr  is  the  Prandtl  number.  Equation  (3)  is 
equation  (21)  of  Liu  [2000],  except  where  4  —  37/  is  a  typo 
for  the  correct  factor  of  3  —  4 77.  Since  the  observed  gravity 
wave  period  is  ~1.5  hr  and  much  shorter  than  the  inertial 
period  of  ~18.5  hr  at  lidar  location,  the  gravity  wave 
dispersion  relation  can  be  simplified  to  [Fritts  and 
Alexander,  2003] 


where  m  is  the  vertical  wave  number  (m  =  27r/Az).  This 
equation  relates  the  vertical  wave  number  to  the  intrinsic 
phase  speed  of  the  wave  and  the  Brunt- Vaisalla  frequency. 

[s]  Since  the  CSU  sodium  lidar  system  can  simultaneously 
and  continuously  observe  the  vertical  profiles  of  temperature, 
zonal  wind  and  meridional  wind  with  high  temporal  and 
vertical  resolutions,  we  are  able  to  obtain  the  time  derivative 
of  the  temperature  and  horizontal  wind,  to  be  used  for  the  left- 
hand  side  of  equations  (1)  and  (3).  The  H,  u,  and  N  in 
equation  (1),  T0  in  equation  (3),  and  dii/dz  in  equation  (2) 
could  be  derived  from  observed  background  atmosphere 
before  the  wave  breaking.  In  fact,  if  £  is  treated  as  an 
adjustable  parameter,  there  are  only  4  unknowns  (c,  k, 
Km,  m )  in  the  4  equations,  (l)-(4).  Therefore,  we  can  then 


infer  the  wave  characteristics,  such  as  vertical  wavelength 
(Az),  horizontal  wavelength  (A/,),  horizontal  phase  speed  (c), 
and  wave  period  (A*/c)  as  well  as  the  eddy  diffusion 
coefficient  (Km). 

[9]  Shown  in  Figure  3  are  the  vertical  profiles  of  heating 
rate,  zonal  wind  acceleration,  and  meridional  wind  acceler¬ 
ation  of  the  mean  flow  at  0900  UT.  These  profiles  were 
obtained  by  first  smoothing  the  vertical  profiles  (Figure  1 
but  with  4  km  vertical  resolution  to  reduce  statistical 
uncertainty)  over  1.5  hours  and  then  removing  the  tidal 
components  of  24  hr,  12  hr,  8  hr,  and  6  hr  periods  prior  to 
the  calculation  of  time  derivative,  with  which  the  rate  of 
changes  between  two  smoothed  profiles  at  45  min  before 
and  after  is  evaluated.  There  exist  positive  peaks  for  both 
zonal  wind  acceleration  and  meridional  wind  acceleration 
near  100  km  with  the  values  of  17  m/s/hr  and  11  m/s/hr 
respectively.  The  negative  peak  in  heating  rate  of  1.5  K/hr  at 
100  km  and  the  positive  peak  of  2.6  K/hr  near  96  km  are 
indicative  of  downward  heat  flux  associated  with  wave 
breaking.  The  negative  peaks  below  95  km  in  the  heating 
rate  and  meridional  wind  acceleration  are  due  to  wave 
oscillations,  as  we  saw  a  clear  oscillation  below  95  km 
with  downward  phase  progression  in  temperature  and 
meridional  wind  profiles  (Figure  1). 

[10]  Since  the  onset  of  wave  breaking  is  near  0900  UT,  we 
use  the  1.5  hr  smoothed  vertical  profiles  at  0815  UT  for  the 
calculation  of  background  atmosphere  parameters  before 
wave  breaking.  The  background  temperature,  zonal  wind, 
meridional  wind,  Brunt- Vaisalla  frequency,  and  mean  wind 
shear  at  the  wave  propagation  direction  near  100  km  are 
188.3  K,  —11.6  m/s,  17.1  m/s,  0.015  s-1,  and  0.55  m/s/km, 
respectively.  Using  equations  ( 1 )— (4),  we  are  able  to  deter¬ 
mine  the  gravity  wave  characteristics.  However,  the  wave 
characteristics  along  with  eddy  diffusion  coefficient  deduced 
will  be  sensitive  to  the  choice  of  the  value  for  £  in 
equation  (3).  By  trial  and  error,  it  is  found  that  with  £  =  4, 
the  apparent  wave  period  detennined  from  equations  (1-4) 
is  85  min,  which  is  close  to  the  1.5  hr  period  revealed  in 
Figure  2.  Near  100  km  at  0900  UT,  this  gravity  wave 
propagates  at  ~30°  north  of  east  as  can  be  determined  by 
the  vector  of  horizontal  wind  acceleration.  This  choice  of  £ 
yields  a  wave  intrinsic  period  of  79  min,  an  apparent  phase 
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Table  1.  Characteristics  of  the  Breaking  Monochromatic  Gravity 
Wave3 


C*(3  —  *n)/v 

3.5 

4.0 

4.5 

Apparent  period,  min 

36 

85 

155 

Apparent  phase  speed,  m/s 

11 

20.5 

30 

Intrinsic  period,  min 

31 

79 

145 

Intrinsic  phase  speed,  m/s 

12.6 

22.1 

31.6 

Vertical  wavelength,  km 

5.3 

9.4 

13.5 

Horizontal  wavelength,  km 

24 

105 

278 

Eddy  diffusion  coefficient,  m2/s 

315 

550 

790 

aHorizontal  wave  propagation  direction:  30°  North  of  East. 


speed  of  20.5  m/s  (22.1  m/s  intrinsic),  vertical  wavelength 
of  9.4  km,  and  horizontal  wavelength  of  105  km.  To 
examine  the  sensitivity  of  the  solutions  to  £,  we  also  solve 
the  equations  with  £  =  3.5  and  4.5  and  the  results  are  given 
in  Table  1.  The  apparent  wave  periods  in  these  two  cases 
were  determined  to  be  36  min  and  155  min  respectively, 
different  from  the  observed  period  of  1.5  hr. 

[it]  For  the  3  choices  of  £  (3.5,  4.0,  and  4.5),  the  deduced 
eddy  viscosities  are,  respectively,  315,  550,  and  790  m2/s. 

In  view  of  its  variability,  these  are  all  consistent  with  the 
average  value  of  320  m2s_l  determined  from  the  momentum 
flux  measurements  by  sodium  lidar  [ Gardner  and  Yang, 
1998].  The  choice  of  £  =  4.0  provides  a  constraint  on 
the  Prandtl  number  and  turbulence  localization  measure: 
3—477  =  4  Pr.  If  we  assume  isotropic  turbulence  (77  =  0), 
then  the  Pr  number  should  be  0.75,  and  if  77  =  0.25,  then 
Pr  =  0.5.  For  the  onset  of  the  gravity  wave  breaking  and 
turbulence,  these  values  of  77  and  Pr  are  plausible.  Using 
the  value  of  £  =  4.0  and  observed  background  temperature 
and  scale  height,  we  can  calculate,  as  an  independent  check, 
the  cooling  rate  due  to  wave  breaking  from  equation  (20) 
of  Liu  [2000],  which  resulted  in  ~1.9  K/hr,  consistent  with 
the  observed  value  at  100  km  in  Figure  3. 

[12]  To  extract  amplitudes  and  phases  of  this  gravity 
wave  packet  with  a  period  of  1.5  hr  from  observation,  we  1 
performed  a  linear  least-square  fit  with  mean  plus  1.5  hr  | 
period  wave  on  the  raw  dataset  between  0600  and  0900  UT, 
yielding  the  best-fit  vertical  profiles  of  amplitudes  and 
phases.  By  visually  examining  the  fitting,  we  are  confident 
about  the  fitting  results.  The  vertical  wavelength  derived 
from  phase  plots  is  about  10-15  km,  which  is  consistent 
with  the  estimation  of  the  middle  column  in  Table  1.  The 
temperature,  zonal  wind,  and  meridional  wind  amplitudes 
were  determined  to  be  ^5  K,  10-15  m/s,  and  10-15  m/s  1 
respectively,  near  100  km.  The  horizontal  wind  amplitude  is  | 
thus  obtained  to  be  ^15-21  m/s,  which  is  close  to  the  wave  | 
intrinsic  phase  speed  of  ~22  m/s.  This  lends  further  support 
for  the  saturation  of  the  1.5  hr  gravity  wave  via  convective 
instability  near  100  km,  since  the  linear  theory  suggests  the 
horizontal  wind  perturbation  is  equal  to  the  wave  intrinsic 
phase  speed  when  the  wave  becomes  convectively  unstable 
[Fritts  and  Alexander,  2003], 

[13]  Shown  in  the  Figure  4  are  the  vertical  profiles  of 
Brunt- Vaisalla  frequency  square  (N2)  and  Richardson  num¬ 
ber  (Ri)  between  0800  and  1300  UT  based  on  the  same 
dataset  shown  in  Figure  1.  These  profiles  are  spaced  at 
15  min  intervals,  whose  horizontal  scale  correspond  to 
5  x  10“4s_2  in  N2  or  3  in  Ri.  The  vertical  dotted  lines 
denote  the  instability  threshold  corresponding  to  0  for  N2 


or  0.25  for  Ri.  Within  the  error  bar  indicated,  the  dynamic 
instability  threshold  of  0.25  may  have  been  reached  at  0815 
UT  near  100  km  and  the  convective  instability  threshold  of 
0  was  then  reached  at  both  0830  and  0845  UT.  When  the 
significant  changes  of  both  temperature  and  horizontal  wind 
were  found  near  0900  UT  (Figures  1  and  3),  however,  the  N2 
and  Ri  near  100  km  are  1.0  x  10-4  s-2  and  0.7  respectively. 
The  1.5  hr  period  gravity  wave  could  start  to  break  after 
0845  UT,  and  then  the  turbulence  mixing  locally  generated 
within  the  breaking  gravity  wave  field  could  bring  the  lapse 
rate  to  adiabatic  value  and  make  the  atmosphere  more  stable. 
The  acceleration  due  to  wave  breaking  will  increase  the  shear 
below  the  wave  breaking  zone  and  in  turn  cause  secondary 
instability  [Liu  et  al.,  1999;  Fritts  and  Alexander,  2003].  At 
0915  UT,  the  N2  was  observed  to  be  more  negative  near 
98  km  and  wind  shear  is  also  dramatically  increased.  After 
9:15  UT,  atmospheric  layers  near  convective  or  dynamic 
instability  developed  below  100  km,  with  lower  boundary 
starting  at  ~98  km  at  9:15  UT  progressing  downward  with  a 
speed  of  ~1  km/hr,  in  agreement  with  tidal  phase  propaga¬ 
tion  during  this  time  period  of  2  to  3  hours.  The  time  series 
in  Figure  4  is  consistent  with  initial  wave  breaking  near 
0900  UT  via  convective  instability,  which  then  trigger  and 
influence  the  subsequent  instability  and  mixing  [Fritts  and 
Alexander,  2003];  this  lasted  until  ^1200  UT  in  a  manner 
similar  to  secondary  instability  under  strong  shear  condition 
investigated  by  Isler  et  al.  [1994]. 

4.  Conclusion 

[14]  On  the  night  of  December  3rd,  2004,  the  Colorado 
State  University  sodium  lidar  system  at  Fort  Collins,  CO 


(a)  Vertical  profiles  of  the  Brunt- Vaisala  frequency  square 


(b)  Vertical  profiles  of  the  Richardson  number 


Figure  4.  Vertical  profiles  of  (a)  Brunt- Vaisalla  frequency 
square  and  (b)  Richardson  number  between  0800  and 
1300  UT  on  the  night  of  December  3rd  (UT  day  338),  2004. 
These  profiles  are  spaced  at  15  min  intervals,  whose 
horizontal  scale  also  corresponds  to  5  x  I0~4s  2  in  Brunt- 
Vaisalla  frequency  square  or  3  in  Richardson  number.  The 
vertical  dotted  lines  denote  the  instability  threshold 
corresponding  to  0  for  Brunt- Vaisalla  frequency  square  or 
0.25  for  Richardson  number. 
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observed  significant  horizontal  wind  acceleration  near  1 00  km 
accompanied  by  a  temperature  cooling  at  the  same  altitude  and 
a  temperature  warming  near  96  km  between  0900-0915  UT. 
This  phenomenon  is  consistent  with  the  modeling  prediction 
for  the  interaction  between  gravity  wave  and  mean  flow. 
Employing  the  gravity  wave  linear  saturation  theory  and 
dispersion  relation,  we  deduced  realistic  values  for  eddy 
viscosity  as  well  as  breaking  wave  characteristics  from  the 
observed  heating  rate,  horizontal  wind  acceleration,  and  back¬ 
ground  atmosphere  without  invoking  modeling  simulation.  By 
approximately  equating  the  calculated  wave  period  to  the 
observed  apparent  period,  we  determined  the  simple  relation 
between  Prandtl  number  and  turbulence  localization  measure 
imposed  by  the  linear  saturation  theory. 

[is]  The  1.5  hr  period  oscillation  at  100  km  is  clearly 
seen  in  the  zonal  wind  component  before  0900  UT,  but  not 
after  0900  UT.  The  Lomb  spectrum  analysis  also  showed  the 
similar  scenario  that  the  power  of  gravity  wave  with  1.5  hr 
period  is  strong  with  above  90%  confidence  level  between 
0500  and  0900  UT,  but  very  weak  between  0900  and 
1300  UT,  indicating  the  wave  saturation  or  breaking.  The 
horizontal  wind  amplitude  deduced  from  observation  is 
~  15-21  m/s,  which  is  close  to  the  wave  intrinsic  phase 
speed  of  ~22  m/s.  The  averaged  vertical  wavelength  deter¬ 
mined  from  vertical  phase  plot  is  consistent  with  the  estima¬ 
tion  from  the  linear  saturation  wave  theory.  Therefore,  it  is 
most  likely  that  the  temperature  warming  and  mean  flow 
acceleration  between  0900-0915  UT  are  due  to  the  breaking 
of  an  upward  propagating  gravity  wave  with  an  apparent 
period  of  ~1.5  hr.  The  instability  analysis  reveals  that  the 
1.5  hr  gravity  wave  could  start  to  break  down  after  0845  UT 
via  convective  instability,  and  then  wave  breaking  may  cause 
the  subsequent  instability  and  mixing  after  0915  UT,  in  a 
manner  indicative  of  secondary  instability  under  strong  shear 
condition. 

[  1 6]  To  our  knowledge,  this  may  very  well  be  the  first 
study,  using  observed  simultaneous  temperature  and  hori¬ 
zontal  wind  profiles  not  only  to  reveal  the  presence  of  wave 
breaking,  but  also  to  estimate  the  characteristics  of  the 
breaking  wave  and  the  coefficient  of  eddy  viscosity.  As 
such,  we  have  demonstrated  the  power  of  simultaneous  and 
continuous  observations  of  temperature  and  horizontal  wind 
with  good  temporal  and  vertical  resolutions,  by  instrument 
such  as  the  Colorado  State  University  sodium  lidar,  for  the 
study  of  wave  breaking  in  the  mesopause  region  of  the 
atmosphere.  Studies  of  this  type  will  be  further  enhanced  by 
accompanying  modeling  on  the  one  hand,  and  by  concur¬ 
rent  imager  observations  on  the  other  to  provide  the 
horizontal  characteristics  of  the  breaking  wave,  checking 
the  horizontal  wavelength  deduced  from  the  linear  satura¬ 
tion  theory. 
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